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We introduce a nonperturbative, first principles numerical approach for solving time-dependent 
problems in quantum field theory, using light-front quantization. As a first application we consider 
QED in a strong background field, and the process of non-linear Compton scattering in which an 
electron is excited by the background and emits a photon. We track the evolution of the quantum 
state as a function of time. Observables, such as the invariant mass of the electron-photon pair, are 
first checked against results from perturbation theory, for suitable parameters. We then proceed to 
a test case in the strong background field regime and discuss the various nonperturbative effects 
revealed by the approach. 
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I. INTRODUCTION 

Treating quantum field theory in the nonperturbative 
regime remains a significant challenge. "Basis Light- 
Front Quantization" (BLFQ) [T], which adopts light- 
front quantization and the Hamiltonian formalism, offers 
a first-principles approach to nonperturbative QFT ^[2]. 
Diagonalization of the full Hamiltonian of the quantum 
field theory yields the physical eigenvalues and eigenvec- 
tors of the mass eigenstates. This approach offers new 
insights into bound state properties and scattering pro- 
cesses |3] as well as opportunities to address many out- 
standing puzzles in nuclear and particle physics (Sj ■ 

The BLFQ approach is real-time (as opposed to 
imaginary-time, as used in lattice-QFT) and therefore 
naturally applicable to time-dependent problems. There 
is currently much interest in gauge theories with an ex- 
plicit time dependence introduced by a background field, 
in particular QED in ultra-intense laser fields [3 [7] and 
QCD in strong magnetic fields [SHIT]. In both cases, the 
greatest interest lies in the case for which the fields are 
strong enough to require a nonperturbative treatment 
and this motivates the approach we present here. 

In this paper we introduce time- dependent Basis Light- 
Front Quantization (tBLFQ), which is an extension of 
BLFQ to time-dependent problems in quantum field the- 
ory. In this approach, BLFQ provides the eigenstates of 
the time-independent part of the Hamiltonian. We then 
solve for the time evolution of a chosen initial state un- 
der the influence of an applied background field, which is 
introduced through explicitly time-dependent interaction 
terms in the Hamiltonian. Although we treat a specific 
application in the present work, the method is more gen- 
erally applicable to time evolution even in the absence of 
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external fields where one is simply following the evolution 
of a chosen non-stationary state of the system. 

In this paper we will apply tBLFQ to "strong field 
QED", in which the background field models the high- 
intensity fields of modern laser systems. Such light 
sources now routinely reach intensities of 10^^ W/cm^, 
and there is ongoing research into using intense lasers to 
investigate previously unmeasured effects such as vacuum 
birefringence |12| and Schwinger pair production [13j . 
Within this research field, the use of large-scale numerical 
codes, based on kinetic models, is becoming increasingly 
popular. Such codes have been used to investigate, for 
example, electromagnetic cascades in laser pulses [T31- 
16 . The two main advantages of such approaches are 
that they are real-time, and that huge numbers of parti- 
cles can be treated via particle-in-cell (PIC) simulations. 
However, there exists no first-principles derivation of the 
required kinetic equations from QED. Consequently, this 
approach is based on a forced welding of classical and 
quantum theories, in which particles and photons are 
treated as classical ballistic objects, and QED cross sec- 
tions are added by hand to model instantaneous colli- 
sions. This leads to problems with double-counting and 
the inclusion of higher-order processes. 

Here, we consider an alternative approach. We restrict 
ourselves to low numbers of particles, but we perform 
a fully quantum and real-time calculation within QED. 
Specifically, we will study "non-linear Compton scatter- 
ing" (nCs), in which an electron is excited by a back- 
ground field and emits a photon [T71 [TS]. This is one 
of the simplest background field processes, as there are 
no thresholds to overcome as in, say, pair creation. We 
note that light-front quantization is the natural setting 
for this investigation [HI |^ , since lasers have inherently 
"light-front" properties: all photons propagate on the 
light-front. 

This paper is organized as follows. We provide the 
background to our approach in Sec. [llj followed by de- 
tails of the BLFQ method in Sec. |III| We then introduce 
tBLFQ in Sec. |IV| and provide illustrative numerical re- 
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suits for our first application to non-linear Compton scat- 
tering in Sec. |V] We present our conclusions and outlook 
in Sec. |VI| The Appendices contain a number of useful 
details. 



The formal solution to Q is 



|V;x+)7=r+exp( -- / Vi]\,P;0)i, (6) 



II. BACKGROUND 

Our approach is based on light-front quantization, and 
on a previously developed method called BLFQ [1] [2]. 
We begin here with a brief review of relevant aspects of 
the light-front formalism O [5T] , and an outline, in terms 
of textbook methods, of the calculation which we wish to 
perform. 

Physical processes in light-front dynamics are de- 
scribed in terms of light- front coordinates {x~^ , a;^. 



"), in which x —a 



plays the role of time. 



Hence, quantization surfaces are null hyperplanes given 
by a;^=constant, and on which initial conditions are spec- 
ified. x~=x'^—x^ is the "longitudinal" direction, and the 
remaining two spatial directions are called "transverse" , 
x-^={x^, x^}. The evolution of quantum states is gov- 
erned as usual by the Schrodigner equation, which in 
light-front quantization takes the form 



^^IV';2;+> = ^p-{x+)\^;x+) , 



(1) 



where | ^; x^ ) is the (Schrodinger picture) state at light- 
front time x^ and is the light-front Hamiltonian. Our 
Hamiltonian will contain two parts; ^qed which is the 
full light-front Hamiltonian of QED, and V which con- 
tains interaction terms introduced by a background field, 
so 



p-{x^ 



P, 



QED 



V(x+) 



(2) 



V contains, in general, an explicit time dependence. It 
is therefore natural to use an interaction picture, but we 
must immediately stress two things: first, we are not us- 
ing the usual "free -I- interacting" split of the Hamiltonian 
and, second, we are not working in perturbation theory. 
Instead, the full QED Hamiltonian -Pqed replaces the 
customary "free" Hamiltonian, and V is naturally the 
interaction term. The interaction picture states are then 
defined by 

|V;x+)/ = ei^«-o^^|V;x+) , (3) 

(since ^qed time-independent), and obey 

i^\t.^^)i = \vi{x+)\ir,x+)i, (4) 

in which V/ , "the interaction Hamiltonian in the interac- 
tion picture" , is 

Vi{x+) = e5^QED^+y(a;+)e-i^^QED^+ . (5) 



where T+ is light-front time ordering. Now let us imagine 
that we could "solve" QED and identify the eigenstates 



and eigenvalues of the theory, 
respectively, so 



Call these 



and P, 



P 



P, 



QED 



\P)^P^\P) 



(7) 



Having these covariant solutions, we would then be inter- 
ested in the transitions between such states introduced 
by the background field interactions contained in V. We 
choose the external field (modeling an intense laser) to 
vanish , ^ = 0, prior to x~^ = 0. At x^ = 0, we expand 
a chosen initial state as a sum over QED eigenstates: 



|V;0)/ = ^|/3)c^(0), 

where c^(0) is the initial data such that 
cMO)ee (/3|^;0)/ . 



(8) 



(9) 



We then expand a solution of the interaction picture state 
at later times. 



I il;;x~' 



]cpix+)\P) 



(10) 



in which the coefficients characterize the nontrivial 
part of the state's time evolution induced by the external 
field. Plugging ( 10 ) into Q yields an equation for the 
c/3: 



.dcp{x^ 



= J2{(3\^Vj{x+)\(3')c0.{x+) 



--: Mpp'{x^)cp'{x+). 



(11) 



(Summation notation in the second line.) This is an in- 
tractable infinite-dimensional system of coupled differ- 
ential equations, and it is at this point that one would 
normally switch to perturbation theory in the interac- 
tion V . However, the background fields we wish to treat 
are strong and therefore not amenable to perturbation 
theory. We therefore write down the formal solution to 
( 11 1, which is, regarding as a column vector and Mppi 
as a matrix, both with infinite dimensions, 



c{x'^) ~ 7+ exp 



M c(0) 



(12) 



In our approach, BLFQ provides finite dimensional ap- 
proximate solutions for the eigenstates |/3). In tBLFQ, 
the time evolution in ( 12 ) is performed numerically, be- 



ginning with the initial vector c(0), to find the vector 
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c(a;+). The coefficients c^(a;+) can then be read off, al- 
lowing one to reconstruct the evolved state itself from the 
overlap 



II ■ 



(13) 



In this way we solve equation (111 with initial condi- 
tions 

Let us compare the above to the usual calculation 
of scattering amplitudes in QED. Such amplitudes are 
based on the split of the QED Hamiltonian into a free 
particle Hamiltonian, Pf^^^,, and an interaction. For the 
application here, this split produces an interaction that 
would be the sum of the QED interaction terms, call them 
Vq, and the additional interaction terms introduced by 
the background, V. 

A scattering calculation would begin with an initial 
state which is a free particle state | i ), prepared at x'^= — 
oo. This state would be evolved through all time using 
the S'-matrix operator |22j . 




FIG. 1. An illustration of non-linear Compton scattering. An 
electron enters a laser field, is accelerated, and emits a photon. 
After emission the electron can be further accelerated until it 
leaves the field. 



S= lim es^froc^^T- e e^s^^f.-?^ . (14 

T-j-oo 

and projected onto a final state | / ) , describing free parti- 
cles at 2;'''= + oo. Thus, one obtains the S'-matrix element 



(15) 



We are also calculating "scattering amplitudes", but 
there are two important differences between our approach 
and that based on the S'-matrix. First, we calculate tran- 
sitions based upon the eigenstate basis of QED (for ex- 
ample physical electrons) rather than between free parti- 
cle states. Second, and related, we calculate finite-time, 
rather than asymptotic, transitions between such states. 
For all times before and after the external field acts on 
our chosen state, we have, in principle, the full quantum 
amplitude expressed as a superposition of physical states 
(mass eigenstates of QED). A specific experimental setup 
will then project this full amplitude onto states to which 
that setup is sensitive. 



A. Application: Nonlinear Compton Scattering 

In this paper we apply tBLFQ to the process of sin- 
gle photon emission from an electron accelerated by a 
background field. Taking the background to model an 
intense laser, this process often goes by the name "non- 
linear Compton scattering" and is well-studied in plane 
wave backgrounds [551 - f2H| . An appropriate experimental 
setup would see the (almost head on) collision of an elec- 
tron with the laser, and the subsequent measurement of 
either the emitted photon [57] or electron [55] spectra. 

We begin with an electron at light-front time x"'"=0 
when it first encounters the laser field. The electron may 
be both accelerated (invariant mass unchanged but 4- 
vector altered) and excited (invariant mass changed) by 



the laser field. Excitation produces electron-photon fi- 
nal states. After time Ax+ the background field switches 
off and no further acceleration or excitation may occur. 
This setup is sketched in Fig. [T] for two of the four di- 
mensions in the problem. The natural question to ask 
is how the quantum states of the electron and (emitted) 
photon fields evolve with light-front time x^, and this 
will indeed be studied below. 

While, in principle, there is nothing to stop us includ- 
ing arbitrarily complex background fields, as a first step 
we consider a simple model. The background is turned 
on only for finite light-front time Aa;^, during which it is 
independent of a;+ but inhomogeneous in , 

eA~(x~) = 2meao cos (Z_x~) (16) 
— nieao [exp {il_x^) + cxp [—il_x~)\ . 

where e is the electron charge and TOg is the electron 
mass. We have written out the exponential form of cosine 
to highlight that the field both "pushes" and "pulls" par- 
ticles in the longitudinal direction. This field has periodic 
structure in the longitudinal direction with frequency l_ 
and dimensionless intensity parameter ag. (aq — 1 cor- 
responds to an intensity of ~ 10^^ W/cm^ at optical fre- 
quency [6 .) It is uniform in the transverse plane, as for 
plane waves, but unli ke p lane waves is longitudinally po- 
larizecQ The profile ([l6]) describes, in the lab frame, a 
beam of finite duration ^/2A.x'^ propagating along the 
x^ direction. Classically, such a field accelerates charges 
in the x~ {x^) direction as time x'^ (a;") evolves. The 



^ It is common to model the focus of a laser field as, for example, a 
time-dependent electric field. Such fields, like ( |16| , do not obey 
Maxwell's equations in vacuum. Pulse profiles which exactly 
obey Maxwell's equations do exist |29| . but they are beyond the 
scope of this first treatment. 
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accelerated charges subsequently radiate, see Fig. [Tj and 
it is the quantum version of this radiation which wc will 
investigate below. 



III. BASIS LIGHT-FRONT QUANTIZATION 
(BLFQ) 

We are interested in how eigenstatcs of the full QED 
Hamiltonian ^^qed evolve due to interactions with a 
background field. (This is analogous to, but clearly not 
the same as, studying transitions between bare states in- 
duced by perturbative QED interactions.) To begin, wc 
must therefore find the eigenstatcs of QED, for which we 
must adopt an approximation. 

The method we use to construct the approxi- 
mate eigenstatcs is Basis Light-front Quantization, or 
BLFQ [2 [2] . This is a numerical method for calculating 
the spectrum of a Hamiltonian, using light-front quanti- 
zation. The idea of finding, for example, the bound state 
spectrum via diagonalization of the Hamiltonian has a 
long history [3 . One well-known approach is discretized 
light-cone quantization [3J [5UH5^ , on which BLFQ is in 
part based. The idea behind BLFQ, and its main advan- 
tage, is that its adopted basis should have the same sym- 
metries as the full QED or QCD Hamiltonian. (BLFQ 
was initially designed for QCD |lj and is supported by 
successful anti-de Sitter QCD methods [33].) This ba- 
sis is therefore not the usual basis of momentum states. 
Usually, the more symmetries the basis captures, the 
less the computational effort that is needed for the solu- 
tions to refiect those symmetries. Because of this, BLFQ 
achieves an accurate representation of the Hamiltonian 
using available computational resources. The construc- 
tion of the BLFQ basis therefore begins with symmetries 
of the light-front Hamiltonian. 

A. Basis construction 

The derivation of the light-front QED Hamiltonian (in 
the presence of background fields) and a list of relevant 
mutually commuting operators, may be found in Ap- 
pendix [A] We do not need the detailed form of these 
operators in order to discuss the three symmetries di- 
rectly encoded in the BLFQ basis. A fourth symmetry, 
transverse boost invariance, (also referred to as trans- 
verse Galilei invariance, [3l[6]) is discussed separately be- 
low as it is not encoded directly in the BLFQ basis but 
is easily accessible with the employed transverse basis. 

The three directly encoded symmetries are 1) Transla- 
tional symmetry in the longitudinal x~ direction. The 
longitudinal momentum operator, P+, therefore com- 
mutes with the Hamiltonian [Pqej^, P^] = 0, and total 
longitudinal momentum is conserved. 2) Rotational sym- 
metry in the transverse plane. This means that the lon- 
gitudinal projection of angular momentum is conserved, 
and the corresponding operator obeys [J^qed' — ^■ 



The operator can be decomposed into two parts for 
each particle species, 

J3 = J3 + J3 , (17) 

in which the subscript "o" refers to the longitudinal pro- 
jection of orbital angular momentum, while subscript "i" 
refers to the longitudinal projection of the spin angular 
momentum. This defines the helicity of a particle in light- 
front dynamics. 3) Charge conservation [-Pqed,Q] = 0, 
where Q is the charge operator with eigenvalue equal to 
the net fermion number Nf. 

The existence of these conserved quantities means that 
the QED eigenspace can be divided up into "segments" , 
which are groups of eigenstatcs with definite eigenvaluef]^ 
P+ cx K, J3 = Mj and Q = Nj. The full spectrum of 
QED is the sum of all such segments. 

The BLFQ basis is chosen to respect these symmetries. 
The essential point is that each basis state, call them | a ) , 
is an eigenstate of the three operators introduced above, 
with the eigenvalues, 

{J^P+,Q}\a) = {M,,K,Nf}\a) . (18) 

Therefore, each state belongs to one and only one seg- 
ment. As a consequence, the BLFQ basis divides into 
segments, and the QED Hamiltonian Pq^j^ is accord- 
ingly block- diagonal in the BLFQ basis. As will be out- 
lined below, this structure allows for a large reduction in 
(numerical) complexity in bound-state calculations. 

The BLFQ basis states are built for each Fock-sector 
(of free-particle states) by allowing the particles to oc- 
cupy orthonormalized modes of a single-particle basis 
that facilitates implementation of the full symmetries. 
The many-particle basis states in each Fock-sector are 
therefore direct products of single particle states, writ- 
ten |Q;),so|a)=®|Q;)in general. It clearly remains to 
specify the details of the single-particle states. 

The single particle basis states are chosen to be two- 
dimensional harmonic oscillator ( "2D-H0" ) states in the 
transverse direction and discretized plane waves in the 
longitudinal direction. This is one choice (among many) 
that facilitates implementation of the symmetries men- 
tioned above. We note in passing the contrast with treat- 
ments of the transverse degrees of freedom in a discretized 
two-dimensional plane wave basis where the orbital pro- 
jection symmetry is lost. We also note the freedom to 
choose another orthonormal basis in the transverse space 
using cylindrical coordinates that may be better for some 
applications. 

Each single particle state carries four quantum num- 
bers, 

a — {k, n, m, A} . (19) 



^ Throughout, K is a natural number. The constant of propor- 
tionaUty is explained below, see equation ||20| 
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The first quantum number, k, labels the particle's longi- 
tudinal momentum. For this degree of freedom we em- 
ploy the usual plane- wave basis states, i.e. eigenstates of 
the fr ee-field longitudinal momentum operator P+, see 
(A21), with corresponding eigenvalues p~^. In this pa- 
per, we compactify x~ to a circle of length 2L. We im- 
pose (anti) periodic boundary conditions on (fermions) 
bosons. As a result, the longitudinal momentum in 
our basis states takes the discrete values 



L 



(20) 



where the dimensionless quantity k~l, 2, 3,... for bosons 
(neglecting the zero mode) and k — ^, |, | for fermions. 
In particular, we have for the laser I* = ^fcias where fcjas 
is a natural number. For convenience, throughout this 
paper we take L = 27r MeV~i so that k can be interpreted 
as the longitudinal momentum in units of MeV. 

The next two quantum numbers, n and m, label the 
degrees of freedom in the transverse directions. As men- 
tioned above we take the transverse components of our 
single particle states to be eigenstates of a 2D-H0 which 
is defined by two parameters, mass M and frequency 
51. (See below for the characteristic scale of the oscil- 
lator, which depends only on a combination of these pa- 
rameters.) These eigenstates are labelled by the quanta 
of the radial excitation, n, and the angular momentum 
quanta, m. The eigenstate carrying these numbers has 
HO eigenenergy 



En,ra = (2n + \m\ + l)n 



(21) 



Since they are not eigenstates of the transverse momen- 
tum operator P-^, the BLFQ basis elements mix states 
with the same intrinsic motion but with different trans- 
verse center-of-mass momenta. This is the price we 
pay for employing the 2D-II0 states as single particle 
basis states in the transverse plane. We may employ, 
when needed, a Lagrange multiplier technique to en- 
force factorization of the transverse center-of-mass com- 
ponent of the amplitude from the internal motion com- 
ponents following techniques used in non-relativistic nu- 
clear physics [331 131] . We may also work with alternative 
coordinates chosen to achieve factorization 

The final quantum number. A, labels the particle's he- 



licity, which is the eigenvalue of Jf , see ( 17 1 . The electron 



(photon) helicity takes values A = ±1/2(A = ±1). 

We present only selected essentials of our method; 
more details of the basis states may be found in Ap- 
pendix |Bj We note here that our transverse modes de- 
pend only on the combination b :— \/Mfl (and not on 
M and Q individually). This is a free parameter which 
must be chosen. Since our goal is to design a basis which 
matches as closely as possible the symmetries of the QED 
Hamiltonian, we note that there is only one mass scale in 
QED, and that is the physical electron mass rUe. A sensi- 
ble choice for our 2D-H0 parameter is therefore]^ 6 — rrie, 



and we adopt this throughout. 

Now, to see why this choice of basis is suited to 
light-front problems, we relate the single particle quan- 
tum numbers {fc, n, m. A} to the segment numbers of the 
states a. So, consider a multi-particle state | a ) — (8)| a ), 
which belongs to a particular segment and is an eigen- 
vector of P"^, J^, and Q with eigenvalues K, Mj and Nf, 
respectively. If fc', m', n-^'', and A' are the quantum num- 
bers for, respectively, the longitudinal momentum, lon- 
gitudinal projection of angular momentum, net fermion 
number and helicity of the Z*^ particle in the state then, 
summing over particles I, we have 



M, = Mt + S 



(22) 

(23) 
(24) 



(The single particle net fermion number is 1 for e, -1 
for e and for 7.) We see that the basis states | a) are 
eigenstates of and Jf individually, with eigenvalues 
Mt and S. Note, though, that it is the sum Mj which is 
conserved by the light-front QED Hamiltonian. 

While each basis state belongs to one and only one 
segment, it is clear that the basis states | a) themselves 
are not eigenstates of QED (written as | /?)). These must 
still be constructed by diagonalizing Pqed ^^^^ basis. 
For example, the physical electron eigenstate | f3 ) = |ephys) 
can be expanded as 



-phys , 



E 



a){a\ 



-phys / 



(25) 



in which both the eigenstate on the left and the basis 
states on the right belong to the same segment. Diag- 
onalizing the Hamiltonian in our basis would yield the 
coefhcients {a\P), and hence the physical states \f3). In 
order to do this, though we need to be able to implement 
our basis numerically, which requires some truncation. 
We turn to this now. 



B. Basis reduction 

Since a quantum field theory contains an infinite num- 
ber of degrees of freedom, reduction of the basis space is 
necessary in order for numerical calculations to be feasi- 
ble. For us, this reduction takes place both in the basis 
states retained (exploiting symmetries) and in the Fock 



In Fock sectors with n particles the effective 2D-HO parame- 



ter for the center-of-mass motion is fej^™ = V nMQ = hy/n, i.e., 
y/n times of that for single-particle states. Thus, in order to 
match the center-of-mass motion across different sectors as re- 
quired by QED vertices, we adopt sector-dependent 2D-HO pa- 
rameters bn = b/y/n for Fock sectors with n particles, where 
b = me is the 2D-HO parameter in the one particle sector. 



6 



space itself (i.e. we retain only certain sectors and imple- 
ment regulators). 

The first type of reduction is called "pruning" , in which 
we exclude basis states which are not needed for desired 
observables. The pruning process is lossless, in that it 
does not lead to loss of accuracy in the desired observ- 
ables. For example, in bound state problems, one is typi- 
cally interested in states with definite Nf and AIj . Com- 
bining this with the longitudinal boost invariance inher- 
ent to light-front dynamics, one can choose K based on 
the desired "resolution" for the longitudinal momentum 
partition among the basis particles Thus, one only 
needs to work in a single segment of the QED eigenspace, 
neglecting the others, without loss of information. From 
here on we write "BLFQ basis" to mean the basis of a 
single segment. 

Pruning alone is not enough to reduce the basis space 
to finite dimension, however, since even a single segment 
contains an infinite number of degrees of freedom. To 
further reduce the basis dimensionality we need to per- 
form basis truncation, which unavoidably causes loss of 
accuracy in calculating observables. Basis truncation is 
implemented at two levels. 

i) Fock- sector truncation. Consider the physical elec- 
tron state. This has components in all Fock-sectors with 
Nf = 1, which we write schematically as 



kphys) = a|e) + b\e-f) + c|e77) + d\eee) 



(26) 



Included in this series are, for example, the bare elec- 
tron I e) and its photon-cloud dressing, | 67), | 677) etc. 
Together, the bare fcrniion and its cloud of virtual parti- 
cles comprise the observable, gauge invariant electron, as 
originally described by Dirac [37H39] . We implement ba- 
sis truncation by assuming that higher Fock-sectors give 
(with an appropriate renormalization procedure imple- 
mented) decreasing contributions for the low-lying eigen- 
states in which we are mostly interested. (One moti- 
vation for this is the success of perturbation theory in 
QED). In this first paper, we make the simplest possi- 
ble nontrivial truncation, which is to truncate our Fock- 
sectors to |e) and \e^). Thus, in this truncated ba- 
sis, the physical electron state would be given by only 
the first two terms of (26). This is enough to calcu- 



late physical wavefunctions accurate up to the first-order 
of the electromagnetic coupling a. Due to its simplicity 
this Fock sector truncation has been typical of light-front 
Hamiltonian approaches such as Refs. 130] though an 
extension to include the 2-photon sector has been success- 
fully implemented in solving for the electron's anomalous 
magnetic moment 41 . 

ii) Truncation within Fock-sectors. Fock-sector trun- 
cation is still not enough to reduce the basis to finite 
dimension; each Fock particle has an infinite number of 
(momentum) degrees of freedom. In BLFQ, truncations 
of the longitudinal and transverse degrees of freedom are 
realized separately, and differently. 

Truncation of the longitudinal basis space is realized 
through the finite size of the x~ direction. By impos- 



ing (anti-) periodic boundary conditions, the longitudi- 
nal momentum k for single particles can only take dis- 
crete values, see (20). Therefore, in a given segment with 



total longitudinal momentum K, only a finite number 
of longitudinal momentum partitions is available for the 
particles in the basis states, since each particle's momen- 
tum must obey Q < k < K and all the fc's must sum to K. 
For segments with larger K, more partitions of longitu- 
dinal momenta among particles are possible, allowing for 
a "finer" description of the longitudinal degrees of free- 
dom. Thus, K also regulates the longitudinal degrees of 
freedom; bases with larger K have simultaneously higher 
ultra-violet (UV) and lower infra-red (IR) cutoffs in the 
longitudinal direction. 

Now consider the transverse part. Recalling from 
above that the transverse states are eigenstates of a 2D- 
HO, with energies (21), we define the total transverse 
quantum number for multi-particle basis states | a ) as. 



I 



2ni + \mi\ 



1 



(27) 



where the sum runs over all particles in the state. This 
number is used as the criterion for transverse basis trun- 
cation; all the retained basis states satisfy 



(28) 



for some chosen A^max- Physically, this simply corre- 
sponds to restricting the total 2D-II0 energy (summed 
over all particles). A'max is specified globally across all 
Fock-sectors to ensure that the transverse motion in dif- 
ferent Fock-sectors is truncated at the same energies. 
As shown in Appendix |B 1[ N^nux determines both the 
UV and IR cutoffs for the transverse basis space, see 

also gmsi. 

This brings us to the end of our discussion on the 
BLFQ basis itself, so let us summarize the approach so 
far. The eigenspace of QED breaks up into segments, 
labelled by K, Mj and Nf. The BLFQ basis is a ba- 
sis of states for such a segment, with each basis element 
carrying the same three quantum numbers as the seg- 
ment itself. The basis elements themselves are collections 
of Fock particle states. For each Fock particle, 2D-II0 
states/plane- waves are employed to represent the trans- 
verse/longitudinal degrees of freedom. The Fock particle 
states carry four quantum numbers, k, n, to. A, see above. 
A complete specification of a BLFQ basis requires 1) the 
segment numbers K, Mj, Nf, 2) the parameters b and L 
pertaining to the transverse oscillator basis and length of 
the longitudinal direction, respectively, and 3) two trun- 
cation parameters, namely the choice of which Fock sec- 
tors to retain, and the transverse truncation parameter 
Nn-Lux- (Recall, K automatically serves as a longitudinal 
truncation parameter because x~ is compact.) 

Such a basis is finite dimensional. It is then a straight- 
forward matter to diagonalize the QED Hamiltonian in 
the BLFQ basis. This yields, as well as the eigenvalues of 
the Hamiltonian, a representation of the physical states 



of QED in terms of the BLFQ basis, as in (25) 
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We end this section with a few words on renormahza- 
tion. Our focus in this paper is on an initial exploration 
of tBLFQ, and we neglect the necessary counter-terms 
when writing down our Hamiltonians. (Hence, we adopt 
physical values for the electron mass and charge.) Renor- 
malization within the BLFQ framework is possible, via a 
sector-dependent scheme [301 HIHiS] . For an application, 
see [15] , in which the scheme is implemented for the QED 
Hamiltonian; the resulting electron anomalous magnetic 
moment agrees with the Schwinger value to within 1%. 



IV. TIME-DEPENDENT BASIS LIGHT-FRONT 
QUANTIZATION (tBLFQ) 

Now that we have the physical states of QED, we turn 
to the transitions between them as caused by an external 
field. Our Hamiltonian now consists of two terms. 



=P, 



QED 



(29) 



in which the new term V comprises the interactions intro- 
duced by the background (laser field), just as in ([2|. See 
Appendix \K\ for the explicit form of the new interactions. 

As discussed above, only a single segment of states is 
needed to address bound-state problems. The presence 
of background field terms V means, in general, that the 
full Hamiltonian will not posses the symmetries asso- 
ciated with conservation of longitudinal momentum [K) 
and longitudinal projection of total angular momentum 
(Mj). {Net fermion number is not affected, of course.) In 
other words, the background field can cause transitions 
between QED eigenstates in different segments. In order 
to account for this, the BLFQ basis must be extended to 
cover several segments. We refer to a collection of multi- 
ple BLFQ basis segments with different iiT's and Af^'s as 
the "extended BLFQ basis" . 

In fact, since our particular choice of background 



field ( 16 1 only adds longitudinal momentum (and light- 



front energy) to the system, the transverse degrees of 
freedom remain untouched, and the symmetry associ- 
ated with holds even with the laser field switched on. 
Therefore, for our current example, we only need to in- 
clude segments with different total longitudinal momenta 
K. 

We therefore begin by applying BLFQ to ^qed ™^ each 
segment, finding the physical states in that segment and 
representing them as in (25 1. The combination of all 
such eigenstates from all the segments forms the "tBLFQ 
basis", which is a basis of physical eigenstates of QED. 
From here on we will write | a ) to represent the extended 
BLFQ basis, and | /3 ) to represent the tBLFQ basis of 
physical states. See Fig. [2] for an illustration of the two 
different bases and the relationship between them. 

In constructing the tBLFQ basis, one needs to specify 
the number of the segments to be included. Larger, less- 
truncated, basis spaces yields more realistic and detailed 
descriptions of the underlying system. The price we pay 



for increased basis dimensionality is of course increased 
computational time. 

We have reached the stage at which we have an (ap- 
propriate) set of physical eigenstates of QED. We now 
describe the preparation of the initial state, and its evo- 
lution in time under the Hamiltonian ( 29 1 . 



A. Initial state preparation 

In perturbation theory, scattering calculations take ini- 
tial states to be eigenstates of the free part of the Hamil- 
tonian, following the usual assumption of asymptotic 
switching, see though [471I3H]- In our calculations, ini- 
tial states are taken to be physical eigenstates of QED. 
For our nCs, process, for example, the initial state is a 
single physical electron with longitudinal momentum Ki. 
This state can be identified as the "ground state" of the 
QED Hamiltonian -Pqeq in the segment Nf=l, Mj = ^ 
and K=Ki (since there is no other state in that segment 
with a lower energy). In the tBLFQ basis, which is just 
the set of eigenvectors of QED, this initial state is triv- 
ially defined. 



B. State evolution 

Recalling the discussion in Section |lTj our initial state 
evolves, in the interaction picture, according to 



IV'; a; 



= T+e 



(30) 



in which \'ip]Q)i is the initial state, equal to a cho- 
sen eigenstate of QED (or a superposition thereof). In 
general, the interaction operator Vj will not commute 
with itself at different times. We decompose the time- 
evolution operator into many small steps in light-front 
time x'^ , introducing the step size Ja;"*", 



;Viixt)6x^ 



(31) 



in which each square bracketed term is a matrix, and 
we let each of these matrices act on the initial state se- 
quentially. Between each matrix multiplication we in- 
sert a (numerically truncated) resolution of the identity, 
so that the evaluation of (31) amounts to the repeated 
computation of the overlaps 

(/?' /3 ) = (/?' |y I /?) exp [| (P^- - P^ )x+], (32) 

in which the are the previously solved eigenergies of 
Pqgj-,, and their presence follows from Eq. (j7|. In order 
to calculate the left hand side of ( [32| in our numerical 
scheme, it is simpler to first calculate the phase factor 
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a ) = 




= 1^ 



-A' etc 



|e> ~ |fc^n^m^A'=)|^,^J 
2) I £7 ) ^ I k\ n^ m^ y)®\k\n\m\V)\ 



lep'^y^; 



I k—K etc 



FIG. 2. The BLFQ and tBLFQ bases. On the left, the extended BLFQ basis | a). This is a collection of bases in different 
segments, each segment labelled by Mj and Nf. (Since nothing in our theory changes net fermion number, all segments 
of interest have fixed Nf = 1, in our case.) The states in each segment are bare states. Two such states, a bare electron and 
a bare electron + a photon, are illustrated. The BLFQ procedure diagonalizes the Hamiltonian in each segment. The basis 
states in I a ) are then rearranged into eigenstates | /3 ) of the QED Hamiltonian, shown on the right. These are the tBLFQ 
basis states. 



and then calculate the remaining overlap in terms of the 
(extended) BLFQ basis, as follows: 



[P'\V\P)^Y.^P'\a'){a'\V\a){a\l3) 



(33) 



The resulting interaction picture matrix elements are the 
elementary building-blocks for evaluating all observables. 

For our particular choice of background field, the struc- 
ture of the matrix elements between BLFQ basis elements 
I a ) is simple. The interaction terms introduced by the 
chosen background do not contain the quantum gauge 
field (see Appendices |A] and |e] for details), and there- 
fore do not directly connect different Fock sectors; ma- 
trix elements of the type (a'(e) \V\ Q;(e7) ) are therefore 
all zero. (Physically, the only direct effect of the cho- 



sen background field is to cither increase or decrease the 
longitudinal momentum k of an electron by /cias-) Ma- 
trix elements between the same Fock sectors (in our case 
{a' {e)\V\a{e)) and {a' {e^)\V\a{e^))), on the other 
hand, are nonzero. If a and a' label two Fock electron 
states, then one finds for example 



(34) 



in which S* is the Kronecker delta. The magnitude of 
the matrix element is proportional to the field intensity 
flQ. It is the sum of two terms, originating in the two 
exponentials in (16 1. Each term is the product of two 
Kronecker deltas. The first delta conserves all quantum 
numbers between the states except for the longitudinal 
momentum (since that is all that our background field 
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alters) . The second delta fixes the difference between the 
k values of the basis elements to be A;' = A: ± fciasor! this 
is simply the "conservation" of longitudinal momentum 
among the initial and final electrons, in that any added 
energy-momentum must come from the laser field. 



a systematic study of nCs using the laser matrix ele- 
ments obtained from BLFQ, in Section VB For inter- 



ested readers we present the full details in the numerical 
calculation for the nCs process (in a "minimal" basis) in 
Appendix. [P] 



Numerical Scheme 



A. Comparison of laser matrix elements 



A direct implementation of Eq. ( |3l[ ) leads to the so- 
called Euler scheme which relates the state at x^+Sx'^ 
to that at x'^; this scheme is however not numerically 
stable (since it is not symmetric in time) and the norm 
of the state vector \ 'ijj;x~^) increases as time evolves, see 
Ref. [in] . We therefore adopt the second order difference 
scheme MSD2 [SD' , which is a symmetrized version of the 
Euler scheme relating the state at x~^+6x'^ to those at 
x'^ and x^—6x^ via 

\i};x++5x+)i 

= I V; X+-SX+ )j + (e-*^^^-+/2 _ e'^^^-+/2)| ^; ^+ ), 



The laser matrix elements 



I x+-6x+ )i ~ iVi{x^)&x'^\ i;; x+ ) 



(35) 



It can be shown that the MSD2 scheme is stable, with the 
norm of the states conserved, provided that Vj-maxSx^ < 
1, where V/;inax is the largest eigenvalue of V/ [35] • This 
requirement imposes an upper limit on the step size Sx~^ . 
Further limits on 6x^ will be discussed below. 

(Note that in order to provide sufficient initial condi- 
tions for the MSD2 scheme, we use the standard Euler 
scheme to evolve the initial state one half-step forward, 
generating | ip; 6x~^ /2 )/. Then we use the MSD2 scheme 
to evolve \'4'',Sx~^/'2)i an additional half-step forward, 
generating | tp; fa"*" )/. With both | ■(/;; 0)/ and | Sx~^ )/ 
available the MSD2 scheme is ready to generate \ip;x~^ ) 
at subsequent times, in time steps of Sx~^ .) 

This concludes our discussion of the principles behind, 
and the method of application, of BLFQ and tBLFQ. 
The reader interested in more details is referred to Ap- 
pendix [C] for the (analytic) representation of states and 
operators in the BLFQ basis, and to Appendix [D] for a 
worked example of the construction of a small, simple 
BLFQ basis, diagonalization of the Hamiltonian and an 
example tBLFQ calculation. 

In the next section we turn to the results of our calcu- 
lation of the nCs process. 



V. NUMERICAL RESULTS 

In this section we present numerical results for non- 
linear Compton scattering (nCs), computed in the 
tBLFQ framework. Since the laser matrix elements 
{P'\V\l3) play an important role in the numerical re- 
sults, we first check them against those from light-front 
perturbation theory, in Section V A We then perform 



[P'\V\f3) 



(36) 



are calculated in the BLFQ framework from the wave- 
functions of I /3 ) and | /?' ) found from diagonalizing 



P, 



QED- 



Due to the small value of the electromagnetic 
coupling a — e'^/{ATr) these wavefunctions can also be 
calculated in perturbation theory, and we will use this to 
check the BLFQ procedure. 

Let us begin with the perturbative calculation of the 
matrix element (36). The background field enters only 
as an operator sandwiched between the states. What 
we must do is to construct the QED eigenstates | /3 ) . 
This can be achieved using ordinary, time-independent 
perturbation theory. To be concrete we will take \ l3) — 
I Gphys ) , the physical electron, and | /?' ) = | e7scat ) 7 the 
electron-photon scattering state. We will work to first 
order in the coupling. So, if | j ) is a complete set of 
eigenstates of the free light-front Hamiltonian Pf7ooi and 
the QED interaction linear in e is Vq (see the first line 
of ( A19| )), then the physical electron can be written, to 
first order, 



I Gphys ) 



{j\VQ\e) 



Similarly, the physical electron-photon state is 

(j|VQ|e7) 



I e7scat ) = I 67 ) - ^ \j\ 



(37) 



(38) 



The matrix element ( 36 ) is therefore approximated in 
perturbation theory by 

(e7scat |V"|ephys) = {e-f\V Vgle) 

-NFrccl^) -'free 

+ {ei\VQ _ ^ ^ —V\e) , 



^f;ee(e7) - Pi 



free 



(39) 



in which we have written a hat over the operator Pf^oe to 
distinguish it from the eigenvalues Pf,-^ci^) and Pf^eeiej). 
Note that |e) and |e7) are eigenstates of Pf7ee' '^^ 
order to compare with the BLFQ calculation we need to 
evaluate the matrix elements (39 1, and hence the states, 
in the BLFQ basis \a). The calculation is uninstructive, 
so we simply present the result in Appendix [E] 
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Now, how do we compare this with a BLFQ calcula- 
tion? We begin by constructing a basis containing only 
two (i4r-)segments, K ~ Ki and K = Ki + fcias, using 
the same parameters as in the perturbative calculation. 
In the K — Ki segment we retain only the single elec- 
tron (ground) state; this acts as the initial state. In the 
K = ifi-|-fcias segment we retain only the electron-photon 
(excited) states. Such a basis, while heavily truncated, 
is all that is required for comparing the results of BLFQ 
and tBLFQ with the perturbative result ( 39 1 . 



In general we would expect that the two matrix el- 
ements will match in the case of small QED coupling 
a. However, in our truncated Fock-space (with only the 
one-electron and one-photon-one-electron sectors) there 
is one further source of potential discrepancy, for the fol- 



lowing reason: the perturbative matrix elements ( 39 1 are 



calculated in the (complete) momentum basis first and 
then projected onto the initial (|e)) and final (|e7)) 
states in the BLFQ basis. The BLFQ matrix elements, 
on the other hand, are calculated in a truncated basis 
space throughout. In the language of perturbation the- 
ory, there exists extra truncation effects in the BLFQ 
matrix elements between the "propagator" and the QED 
vertices, cf. Eq. (39). Due to this "intermediate" basis 



truncation, exact agreement can only be expected in the 
continuum limit (iVmax — > oo). 

Direct comparison as a function of A'max is difficult, be- 
cause as iVmax increases the spectrum of ^qed changes 
(more states appear in the spectrum) and it becomes dif- 
ficult to keep track of the A^max dependence for specific 
matrix elements. We will now look at a test case, the nCs 
process with perturbative and nonperturbative matrix el- 
ements as inputs, and compare their predictions for the 
population of various tBLFQ basis states as a function 

of A'^max- 

We take Ki=1.5, A^max^lG, and parameters ao=10, 
fcias=2 for the laser profile ( 16 1. Thus the laser can cause 



transitions between the K—1.5 and K=3.5 segments. In 
general, we expect the truncation error between pertur- 
bative and nonperburbative matrix elements to dimin- 
ish as iVmax increases and more basis states are used. 
However, since we are neglecting various renormalization 
counter-terms, at sufficiently large A'max, high order and 
divergent loop effects may lead to further discrepancies 
between the (leading order) perturbative and nonperbur- 
bative matrix elements. Contributions from these loop 
effects are proportional to (higher-than-leading) powers 
of a and generally increase with the ultraviolet and/or 
infrared cutoff iVmax- Since at this stage we are inter- 
ested only in verifying the decreasing truncation error 
as iVmax increases, we suppress here the contribution of 
loop effects by artificially reducing the coupling constant 
so that eV(47r) ^1/13700. 

At a;+=0 we switch on the laser field and evolve the ini- 
tial single electron state according to Eq. ( [SS] ) . The popu- 
lation of tBLFQ basis states (the probabilities |c^(a:'^)p) 
in the K=3.5 segment, as a function of light-front time, 
are shown in Fig. [3j along with the corresponding pertur- 
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FIG. 3. (Color online) Time evolution of the electron sys- 
tem in the laser field (at A'^,nax=16). Upper, middle and lower 
panels correspond to exposure time a;+=0, 100, 200 MeV~^ 
respectively (the laser field is switched on at x'^—0). Each 
dot on these plots corresponds to a tBLFQ basis state | /3 ) in 
K=3.5 segment. Y-axis is the probabifity, jc^j^, for each ba- 
sis state and x-axis is the corresponding invariant mass, Mp. 
Green (red) dots are results based on laser matrix elements 
evaluated nonperturbatively (perturbatively). Note that the 
electromagnetic coupling constant Q=e^/(47r) is reduced to 
1/13700, see text for details. 



bative results. Different tBLFQ basis states are distin- 
guished by their respective invariant masses, Mp, as de- 
fined in Eq. ( D2 ) . As time evolves the probability for the 

-0.511MeV) 



single electron state (with invariant mass 
drops and various electron-photon states (with invariant 
mass above 0.6MeV) in the K=3.5 segment are gradually 
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FIG. 4. (Color online) "Snapshots" of the system at a;"*" — 4.5 
MeV~^ in bases of A^max^S (upper panel), 16 (middle panel) 
and 24 (lower panel). Each dot on these plots corresponds 
to a tBLFQ basis state | /3) in K=3.5 segment. Y-axis is the 
probability, \c(i\^, (on a greatly expanded scale compared to 
Fig. [sj for each basis state and x-axis is the corresponding 
invariant mass, Mp. Green (red) dots are results based on 
laser matrix elements evaluated nonperturbatively (pertur- 
batively). Note that the electromagnetic coupling constant 
a—e^/l-iTr) is reduced to 1/13700, see text for details. 
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FIG. 5. (Color online) Time evolution of the average in- 
variant mass of the electron system. Up, middle and lower 
panels are calculated in tBLFQ basis space with A'^max=8, 16, 
24 respectively. Y-axis is the difference between the average 
invariant mass (M) of the system at and that of a sin- 
gle electron rrie. X-axis is the (lightfront) exposure time . 
Green (red) dots are results based on (non)perturbative laser 
matrix elements. Note that the electromagnetic coupling con- 
stant a—e^ /{An) is reduced to 1/13700, see text for details. 



populated. At x+=100 and 200 MeV ^, a peak structure 
is seen around the invariant mass of 0.74 MeV. This can 



be understood as foUovifs: because our laser profile ( 16 ) 



is only trivially dependent on light-front time, in that it 
switches on and off but is otherwise constant, in the infi- 
nite time limit only transitions between basis states with 
the same light-front energy can accumulate (the tran- 
sition amplitudes between states with unequal energies 



oscillate with a period inversely proportional to their en- 
ergy difference) . In this case the light-front energy of the 
initial (single electron) state is P~=ml/ Ki=0. 17 MeV, 
basis states in the K=3.5 segment with invariant mass 
around Mpk=y^ P' 17 x 3.5=0.78 MeV will thus 
accumulate and form a peak. The full peak develops over 
longer times and is located at approximately 0.8 MeV, 
independent of A'max- 

In order to study the convergence between the pertur- 
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bative and nonperturbative laser matrix elements as a 
function of A^max, we consider snapshots of the system 
at a fixed exposure time (a;+=4.5 MeV~^) calculated in 
bases with increasing iV^ax in Fig. |4j As expected the 
overall agreement between the results from perturbative 
and nonperturbative laser matrix elements indeed im- 
proves systematically as the basis A^max increases. 

As a measure of the energy transfer between the sys- 
tem and the laser field, we calculate the evolution of the 
average invariant mass (M(a;+))=X] of the 
system as a function of the exposure time x"*". The nu- 
merical results calculated in bases with A'max=8, 16,24 are 
compared in Fig. [sj 

As the exposure time increases, the laser field pumps 
energy into the system, and the invariant mass of the 
system increases accordingly, as seen from Fig. [Sj Again, 
as expected the agreement between the results from per- 
turbative and nonperturbative laser matrix elements im- 
proves as the basis size (A^max) increases. 

In the next subsection we will study the nCs process 
systematically in a larger basis space using laser matrix 
elements from the BLFQ approach. 



B. Numerical results for nCs 

With the laser matrix elements checked, we now turn 
to nCs in a larger basis. This basis consists of three 
segments with K={Ki, Ki+kii^, Ki+2k\iis} ■ In each seg- 
ment we retain both the single electron (ground) and 
electron-photon (excited) state(s). The initial state for 
the nCs process is a single (ground state) electron in the 
K—Ki segment. This basis allows for the ground state 
to be excited twice by the background (from the segment 
with K=Ki through to segment with Ki+2k\a.s)- In this 
calculation, we take Ki=l.b and A^max=8, with ao=2 and 
fcias=2. We present the evolution of the electron system 
in Fig. [6j at increasing (top to bottom) lightfront time. 
The initial system is shown in the top panel of Fig. [6j the 
only populated basis state is the single electron (ground) 
state in the K—l.h segment. As time evolves, the back- 
ground causes transitions from the ground state to states 
in the K='i.h segment. Both the single electron state and 
electron-photon states are populated; the former repre- 
sent the acceleration of the electron by the background, 
while the later represent the process of radiation. At 
times a;+=0.2 MeV~^, the single electron stattj^in K—i.b 



* Because we neglect counter-terms, the single electron ground 
states in the K=l.b, 3.5 and 5.5 segments receive increasing 
(negative) mass corrections from loop effects. K works as an ul- 
traviolet and infrar ed reg ulator in the longitudinal direction (see 
discussion in Sect. |III B I and as a result, the calculated value 
for the invariant mass of the K='i.b and X=5.5 single electron 
states is slightly lower than that for _f<'=1.5. In order to prevent 
the invariant mass of the whole system being affected by this 
artifact, we manually set the invariant mass for each iC-segment 
single electron state to the physical mass me. 



becomes populated while the probability for finding the 
initial state begins to drop. In the right hand panel, the 
populated electron-photon states begin forming a peak 
structure. The location of the peak is around the in- 
variant mass of 0.8 MeV, roughly consi stent with the 
expected value of Mp]^i = ^J P- [K- + fciaT)^0.78 MeV, cf. 
the discussion in Section fVAl 

Once the basis states in X=3.5 become populated, 
"second" transitions to the K—b.b segment become pos- 
sible. This can be seen in the third row of Fig. [6j at 
x^—QA MeV^^. In the left hand panel, one sees that 
the probability of the electron to remain in its ground 
state {K=l.^) is further decreased, the probability of it 
being accelerated (to K='i.b) is increased, and that the 
K—b.h single electron state becomes populated. In the 
right hand panel, the electron-photon states in the i4r=5.5 
segment also become populated as a result of the sec- 
ond transitions. A second peak arises here at the invari- 
ant mass of around Afpk2=-\/-P~(i^i + 2fcias)~1.0 MeV 
(distinct from at that ~0.8 MeV, above, formed by the 
K—i.h electron-photon states from the first transitions). 
The peak in the K—b.h segment is at a larger invariant 
mass than that in the K=2).b segment simply because 
the basis states in the if=5.5 segment follow from the 
initial state being excited twice by the background field, 
and thus receive more energy than states in the -fC=3.5 
segment. 

As time evolves further, the probability of finding a 
K=b.h single electron exceeds that of finding a K=?>.b 
segment electron, see the bottom left panels in Fig. |6] 
At this time, a;"*"=0.6 MeV~^, the system is most likely 
be found in the K—^.^ single electron state, with prob- 
ability ^0.6. The probability for finding the 7^=3.5 sin- 
gle electron state is around 0.35 and the initial K=l.b 
electron state is almost completely depleted. In the 
right hand panel, we see that the probability for find- 
ing K—h.b electron-photon states increases with time. 
One also notices that at later times, the probability for 
i4r=3.5 electron-photon states also begin to drop. This 
is because (hke the if =3. 5 single electron) the K=2>.b 
electron-photon states are coupled to the if=5.5 single 
electron; as the probability of the K=h.?) single electron 
state increases, it "absorbs" both the single electron and 
the electron-photon states in the K=i.h segment. At 
x+=0.6MeV^^ we terminate the evolution process, as 
the system is already dominated by the single electron 
state in the maximum ii'-segment. Further evolution 
without artifacts would require bases with segments of 
i5:={7.5,9.5...}. 

This calculation, although performed in a basis of lim- 
ited size, illustrates the basic elements of the tBLFQ 
framework. The acceleration of the single electron state 
and the radiation of a photon are treated coherently 
within the same Hilbert space. 

Since the states | /3 ) encode all the information of the 
system, they can be employed to construct other observ- 
ables. As an example, in Fig. [7] we present the evolu- 
tion of the average invariant mass {M) of the system 
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FIG. 6. (Color online) Time evolution of the single electron system in the laser field. From top to bottom, the panels in 
each row successively correspond to lightfront-time x'^—0, 0.2, 0.4, 0.6MeV~^ (the laser field is switched on at a:^=0). Each 
dot on these plots stands for a tBLFQ basis state. Y-axis is the probability for the tBLFQ basis state |c/3(a;''')p and x-axis 
is its corresponding invariant mass Mp. The panels on the left (with y-axis up to 1.1) illustrate the evolution of the single 
electron (ground) states in K=1.5, 3.5, 5.5 segments respectively and the panels on the right with y-axis "zoomed-in" show 
the evolution of various electron-photon (excited) states. The electromagnetic coupling constant a—e^/l^Tr) is 1/137. 



as a function of time. The increase of the invariant mass with time reflects the fact that energy is pumped into the 
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FIG. 7. (Color online) Time evolution of the average invariant 
mass of the electron system calculated in tBLFQ basis space 
with A'^niax^S. Y-axis is the difference between the average 
invariant mass (M) of the system at x'^ and that of a single 
electron me. X-axis is the (lightfront) exposure time x"*". The 
electromagnetic coupling constant a=e^/(47r) is 1/137. 



electron-photon system by the laser field. This invariant 
mass can be accessed experimentally by measuring the 
momenta of both the final electron, p^, and photon, pi^ 
in an nCs experiment. The invariant mass can be com- 
pared with the expectation value of (p^ +P^)'^ measured 
over many repetitions of the nCs experiment. Work in 
deriving other observables, such as the cross sections for 
specific electron-photon final states, is in progress. 

In this section we have demonstrated a) the general 
procedure for treating processes nonperturbatively in 
tBLFQ, and b) the accessibility of the full configuration 
(wavefunction) of the system at finite time. 



VI. CONCLUSIONS AND OUTLOOK 



process in which an electron absorbs energy-momentum 
from an intense background laser field, and emits a sin- 
gle photon. In contrast to current numerical approaches 
to strong laser physics, tBLFQ is fully quantum mechan- 
ical and allows us to sec both the acceleration of the 
electron by the background and the creation of a pho- 
ton, in real-time. Note that tBLFQ is also applicable to 
problems without external fields but in which nontrivial 
time-dependence arises from using an initial state which 
is a non-stationary superposition of mass eigenstates. 

Future developments will be made in two directions. 
First, further improvement of tBLFQ itself. The ini- 
tial step is to implement renormalization so that the 
BLFQ representation of the physical eigenspectrum of 
QED can be improved (and then used in tBLFQ cal- 
culations). Currently we are working on implementing 
a sector-dependent renormahzation scheme within the 
BLFQ framework. The inclusion of higher Fock sectors 
in our calculation is also important, as it will not only 
result in more realistic representations of quantum states 
but will also allow for the description of a larger variety 
of processes, e.g., multi-photon emissions. 

The second direction to be pursued is the extension 
of tBLFQ 's range of applications. In the field of intense 
laser physics, the inclusion of transverse (x^), longitu- 
dinal {x~) and time {x'^) dependent structures to the 
background field will be used to more realistically model 
the focussed beams of next-generation laser facilities 29 . 
In addition to intense laser physics, we will also ap- 
ply tBLFQ to relativistic heavy-ion physics, specifically 
the study of particle production in the strong (color)- 
electromagnetic fields of two colliding nuclei. Ultimately, 
the goal is to use tBLFQ to address strong scattering 
problems with hadrons in the initial and/or final states. 
As supercomputing technology continues to evolve, we 
envision that tBLFQ will become a powerful tool for ex- 
ploring QCD dynamics. 



In this paper we constructed a nonperturbative frame- 
work for time-dependent problems in quantum field the- 
ory, referred to as time-dependent BLFQ (tBLFQ). This 
framework is based on the previously developed Basis 
Light-front Quantization (BLFQ) and adopts the light- 
front Hamiltonian formalism. Given the Hamiltonian 
and the initial configuration of a quantum field system as 
input, the system's subsequent evolution is evaluated by 
solving the Schrodinger equation of light-front dynam- 
ics. The eigenstates of the time-independent part of the 
Hamiltonian, found by the BLFQ approach, provide the 
basis for the time-evolution process. Basis truncation 
and time-step discretization are the only approximations 
in this fully nonperturbative approach. One feature of 
the tBLFQ framework is that the complete wavefunction 
of the quantum field system is accessible at any interme- 
diate time during the evolution, which provides conve- 
nience for detailed studies of time-dependent processes. 

As an initial application we have applied this frame- 
work to an external field problem. We have studied the 



ACKNOWLEDGMENTS 

We acknowledge valuable discussions with K. Tuchin, 
H. Honkanen, S. J. Brodsky, P. Hoyer, P. Wiecki and Y. 
Li. This work was supported in part by the Department 
of Energy under Grant Nos. DE-FG02-87ER40371 and 
DESC0008485 (SciDAC-3/NUCLEI) and by the National 
Science Foundation under Grant No. PHY-0904782. 
A. I. is supported by the Swedish Research Council, 
contract 2011-4221. Diagrams created using JaxoDraw 

EH Eg. 



Appendix A: The light-front QED Hamiltonian 

In this section we follow the derivation of the Hamilto- 
nian in 3 , but with an additional background field. The 
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Lagrangian is 



C = --F^,F^'' + ^{ij^D^-m„)^ , (Al) 



is the sum 



in which Z?^ = + ieC^ and — Af^ 
of the background and quantum gauge fields respectively. 
Note that i^^jy is calculated from alone, i.e. there is 
no kinetic term for the background. The equations of 
motion for the fields are 



which defines the current j'^ , and 



(A2) 



(A3) 



The background field appears in the equations of motion 
for the fermion, but not for the gauge field. We now 
analyze these equations in light-front coordinates (x* = 
± x^, and x"^ — 2xzp). We work in light- front ga uge , 
so that A+ = A+ = 0. The = + component of ( [aS] ) 
does not contain time derivatives, and can be written 



-A- 



d^A^ 
d+ 



(A4) 



This is a constraint equation which relates the (non- 
dynamical) field A" to the transverse components A^ 
and the fermion current. Similarly, if we multiply (A3) 



by 7''' on the left, we find a constraint equation for the 
fermion field. Defining first the orthogonal field compo- 
nents 

= \i+r^ , = i7"7+* , (A5) 

the constraint equation may be written 
1 



2id- 



(A6) 



Hence, the field is non-dynamical and can be ex- 
pressed in terms of the dynamical field 5'+. We now turn 
to the construction of the Hamiltonian. The conjugate 
momenta are 



dC 



dC 

dd^A,. 



(A7) 



and the Hamiltonian P = 2P, is then 



P=j '^^x^dx- Ff'+d+A^ + i^-f+d^^ - C 

= J d'^x^dx- F^'+d+A^ + ^F'^'Ff,^ + i^-f+d+^ , 



in which the first line is the standard Legendre transfor- 
mation, and in the second line we have used the equations 
of motion. It is convenient to add a total derivative to 



the Hamiltonian [3], the term — 9^(F^+A+), and again 
use the equations of motion to write 



(A9) 



In order to complete the transition to the Hamiltonian 
picture we need to eliminate the light-front time deriva- 
tives of the fields in favour of the fields themselves, and 
their momenta. The gauge field terms are simplest. Let 
j, j be transverse indices and define 



{A\A-,A^}:^{0X—,A^} 



(AlO) 



The first line of ( A9 1 then becomes 
1 „ 1 



':F'^F,, 



-F+-F, 



(All) 



using the constraint (A4). The field is that which 
survives the limit e — !■ 0, and is therefore referred to as 



a "free field" . Turning now to the spinor terms in ( A9 1 , 
we have 



(A12) 



and the spinor equations of motion ( A3 ) then give 
1, 



[me - z7^£>-^]7"\['_ 



me 



(A13) 



The first line follows from 7 = 0, we used (A6) in 



the second line and in the third line we commuted 7 to 
the right. In analogy to A, we introduce ^, defined by 



1 

2id 



-[me-i7^5^]7+*+ . (A14) 



Again, this is the field which survives the e — > limit. 
Our final task is to insert (A14) into (A13) and rewrite 



this in terms of only the "tilde" variables. First, the 



(A8) C-free terms of of ( |A12[ ) are: 



za+ 2 io+ 

(A15) 
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Next, we have terms in (A12) which are hnear in C: 



Now, we sum (A13) (A15), (A16I and (A17) to obtain 



the fuU Hamiltonian: we drop the "tilde" on all variables 
from now on, so that one must remember that 



(A16) 



using (A14| in the second line. Note the tilde on j in 



the third line. Finally, we have the terms quadratic in C, 
which are 

-e^^4ej^C^]^[e^^C^]^^ . (A17) 



A- = 2 



d+ 



1 

2id+ 
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(A18) 

Since we are interested in the new interactions introduced 
by the background field, we will separate these out ex- 
plicitly, expanding C ~¥ A + A. (Recall, A has a tilde 
now.) Finally, the full Hamiltonian is 



P- 



i2 i 1 1 7 +"^6 

2 ' 



1 j+ + ^^^^^A^^rA,^ 
2 10+ 



(A19) 



The first line is the QED light-front Hamiltonian, ^^qed- 
The second line contains the new terms generated by the 
background field. We label the terms in Pq^d as T/, T-y, 
Wi . . . W3 respectively. Tf and are the kinetic energy 
terms for the fermion and gauge field respectively. Wi 
is called the vertex interaction, which is responsible for 
photon emission and electron-positron pair-production 
processes. W2 is the instantaneous-photon interaction 
and W3 is the instantaneous-fermion interaction. The 
instantaneous-photon interaction is the light-front ana- 
logue of the Coulomb energy, and its origin is Gauss's 
law (A4) [Sj. The instantaneous-fermion interaction is 



(explicitly) present exclusively in light-front dynamics. 

The additional terms introduced by the background 
field, in the second line of (A19), are the three-point 
background vertex interaction, the instantaneous 2- 
background, 2-fermion vertex, and two, instantaneous, 
1-background, 1-photon, 2-fermion vertices. For the field 



(16 1, the Hamiltonian (A19I contains only a single term 



beyond the ordinary QED Hamiltonian, this term being 



j+A+ = *7+^'^^ 



2¥_^-^+A, 



(A20) 



As the main goal in this work is to introduce the general 
framework of BLFQ (rather than to present new results 
through precise numerical calculations, for which see fu- 
ture articles) we work for convenience with a truncated 
QED Hamiltonian, dropping the instantaneous interac- 
tion terms W2 and W3, proportional to e^. The remain- 
ing Hamiltonian is sufficient for calculating eigenstates 
and eigenvalues to first order in a. 



1. Symmetries of Pqed 

The BLFQ basis explicitly carries three of the sym- 
metries of the QED light-front Hamiltonian -Pqed- ^ 
fourth symmetry, boost invariance in the transverse di- 
rection, is not encoded directly in the basis. However, as 
discussed in the text, due to the choice of HO basis and 
the A^max truncation method, this symmetry is recovered 
by a factorization of the resulting amplitudes into a com- 
ponent for center-of-mass motion times the components 
for internal motion j46j . The three symmetries encoded 
directly in the basis and their operators, which commute 
with -Pqed' fisted below. 

The longitudinal momentum operator is. 



P+ = 



1 



2*^^0+1'+ 



d+A^d+A^ , (A21) 



with j G {1,2}. This commutes with -PqeD' ^^'^ ^'^ over- 
all longitudinal momentum is conserved. 

The longitudinal projection of angular momentum is 
also conserved. The corresponding operator is written 
J^. This can be decomposed into the following four parts. 



j: 



(A22) 



in which the subscript "o" refers to the longitudinal pro- 
jection of orbital angular momentum, while subscript "j" 
refers to the longitudinal projection of the spin angular 
momentum, i.e. helicity. The subscripts / and 7 refer to 
the fermion and photon, as above. In terms of the fields. 
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these four operators are 



73 



73 



7;o ~ 



r - 

7;2 ^ 



i Jdx-d^x^ x'[d+A'd^A' +d+A^d^A^ 



-x^[d+A'd'A' +d+A^d'A^] , 



(A23) 



where S'^ = 



cr-^ 
ct3 



Finally, net fermion number is also conserved. The 
corresponding operator is 



Nf = / dx-d^x-^-^l^^ 



(A24) 



Appendix B: BLFQ harmonic oscillator basis 

Our BLFQ basis elements differ from the usual basis 
of momentum states only in the transverse degrees of 
freedom. In this appendix, we describe the transverse 
structure of our basis elements. 

For the transverse part, the basis elements are eigen- 
states of the following two-dimensional harmonic oscilla- 
tor (2D-H0) Hamiltonian 



TT2d _ 

^HO - 



9 

PI 
2Af 



(Bl) 



in which M and i7 are the mass and frequency of the 
oscillator. The choice of these free parameters will be 
discussed shortly. The characteristic scale of the 2D-H0 
is & = V-A/il will be called the 2D-H0 parameter. The 



eigenstates of (Bl) are labelled by two quantum num- 



bers, n, the principle quantum number characterizing 
the quanta of the radial excitation, and m, the angular 
quantum number characterizing the angular momentum. 
These eigenstates, which we write | nm ) have eigenen- 
ergy En^m = (2n + \m\ + 1)^1. In coordinate space, the 
corresponding wavefunctions can be factorized into a con- 
ventional angular part Xm{4')i 



Xm{(f>) 



/27r 



(B2) 



and a radial part fmnip), as so 



= {x^\nm) = (-l)"zl'"l/„™(p)x„.(0) , 

(B3) 

where {p, (j)) are polar coordinates in the transverse plane, 
= pcoscj) and = psin0. Explicitly, the radial part 



fnmip) is given in terms of generalized Laguerre polyno- 
mials, Llr'(fe^p^), by 



{n + |to|)! 



2 „2\ 



(B4) 



The 2D-II0 wavefunctions iPi ^) satisfies the follow- 
ing orthonormalization condition. 



[ nm I n'm!) 



00 /'27r 



dppd0 



Jo 
= (5"' (5"' 



(B5) 



In general, the 2D-H0 wavefunction ^^^^ (p, (p) is a highly 
oscillatory function with respect to both p and (j). In the 
radial direction, though, the oscillations terminate with 
a steep fall-off to zero at around bp ^ 2^2n -|- |m|. 

One property of these wavefunctions is that their coor- 
dinate and momentum space expressions are very similar. 
To see this, we Fourier-transform to obtain the 

momentum space wavefunction ^J^m(p^), 

= (2i:)~t^{v)xM , (B6) 

in which 

V2 I 



fnrniP) = 



and 



n 

{n + \m\)\ 



im4> 



l™l r\rn\ IP 



/2n 



(B7) 



(B8) 



The coordinate and momentum space wavefunctions (B3 ) 



and ( |B6[ ) differ only in an overall coefficient and in that 
the 2D-H0 parameter b appears in numerators or denom- 
inators, respectively. Note in particular that the wave- 
functions depend only on 6, not on M and f2 individually. 



1. Basis truncation: IR and UV cutoffs 

As discussed in the paper, the BLFQ basis must be 
truncated in order for numerical calculations to be fea- 
sible. One of the conditions used for obtaining a finite 
dimensional basis space is that the transverse degrees of 
freedom of the multi-particle states obeys 



\mi\ 



1 < N„ 



(B9) 



particles 



This restriction also imposes both IR and UV cutoffs into 
our theory, as can be seen from the behavior of the 2D- 
HO wavefunctions. 
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^ 2by/2n + \m\. The 
can be supported by 
irefore around 
which is an ultraviolet cutoff in momen- 



The momentum space HO wavefunction ( B6 1 exhibits 
a sharp faii-ofF at around p± 
maximal transverse momentum 

the basis spaces truncated at N^^^ is therefore around 
pI^"" cx by/W, 
turn space. 

Since the coordinate space wavefunctions are so closely 
related to those in momentum space, we see immediately 
that the same basis states in coordinate space have sup- 
port up to x™'^'^ oc -y/TVinax/^- This translates into an in- 
frared cutoff in the momentum space as p™'" = 1 /x™^'^ oc 
fe/V-^max- The above UV and IR cutoffs are analogous 
to cutoffs of the 3D-H0 that have recently been analyzed 
in low-energy nuclear physics applications [HI |33] . 



Appendix C: Light- front QED in the BLFQ basis 
1. Mode expansion 



in which A* = 7^7^/4 are the usual orthogonal, light- 
front projectors and e{x) is the sign function. 

2. Spin and polarization 



We use the following (standard) Dirac spinors, with 
hehcity A = ±1/2 =t; 



vipA) 



I I \ 



n+ 



P 

p+ 



u{p,i) 



vip,l) 



1 ^ 



\ 



J 

^ 



p+ , 



(C5) 



The mode expansion for field operators in the BLFQ 
basis is, 

,(p^)i;(p,A)e^P-] , (CI) 
■[aa^nm{p-^)et.ip, X)e~^''''' 



a nm \ 



dV 



(C2) 



where p • x = ijp^x — p-^ ■ is the 3-product for the 
spatial components of and x^, and see Eq. (20) for the 



values of p+, which depend on the (anti-)periodic bound- 
ary conditions for (fermions) gauge bosons. The cre- 
ation operators &a, d\ and a\ create electrons, positrons 
and photons (respectively) with quantum numbers a = 
{fc, n, m. A}. They obey the (anti-)commutation relations 



[aa,at,] = {ba,b\,} = {da,d\,} = 6aa' 



(C3) 



With this, and using the explicit forms of the spinors and 
polarization vectors, given below, one can verify that the 
fields obey the standard equal-light-front-time commuta- 
tion relations, 



[Mx),A,{y)\^, 



(C4) 



We use a circularly polarized basis of polarization vectors 
for the photon, with A = ±1 =ti, 



e'^(fc,A)= 0,e^(A), 



2e-L(A) • k 
k+ 



(C6) 



in which the transversal polarization vectors are 
= ^(l,z) and e^(-l) = ^(l,-«). The vec- 
tors are normalized according to 



e''(fc,A)e:(fc,A') = -<5AA' 



(C7) 



The Hamiltonian 



We have now written our free field operators in terms 
of creation and annihilation operators. The next step is 
therefore to express the Hamiltonian in terms of the same 
basis. Hence, we take the operators iCll and insert them 



into (A19). The calculation is lengthy and unenlightcn- 



ing, so we will simply give some example terms. 

First, the kinetic energy term for the fermions, which 
is the first term of (A19) : 
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Tf 



Ax' 



2 ' id+ 



-{b\,ba+dl,de,)x 



y.([ml + (2n + |m| + l)^^]^^ - ^^[^(n + l)(n + H + l^^ ^ + ^n{n + \m\)5i+^\5t 5^ 5^x ) ■ (C8) 



One can similarly obtain expressions for the interaction 



terms in (CI ) 
example 



The vertex interaction Wi becomes, for 



V / * aia2Ct3 \/ P2 

+ ^nrmi{pi)^n2ni2ipi)^nsrn3iP3)u{P3, hh^ep.iP2, \2)uipi, Xl)6''^'^ {pi + P2 ^ PsWa^aa^bsi 

+ six similar terms.] (C9) 



Here ai, 012, 0:3 are the quantum numbers associated with 
the field operators ^, and ^ respectively. The 
functions should be understood as the Dirac delta func- 
tion for the transverse momentum (p^) and the Kro- 
necker delta for the discretized longitudinal momentum 
(p+). The two terms given above cause transitions be- 
tween the I e ) and | 67 ) Fock sectors. The six terms 
which we have not written explicitly do not contribute 
to the calculations in this paper, as they describe tran- 
sitions between Fock sectors which are not present in 
our truncated basis space. The spinor-polarization vec- 
tor contraction part u(jp^, ^3)7^6* (p2, A2)u(pi, Ai) in the 
first term is summarized in Table |l] for different helicity 
configurations. Taking complex conjugates and changing 
labels gives the results with instead for e*. Integra- 
tion over the product of three, highly oscillatory, 2D-H0 
wavef unctions, as in (C9|, would pose a challenge for nu- 
merical calculations. Fortunately, this type of integral 
can be performed analytically by applying the Talmi- 
Moshinsky transformation to the 2D-H0 wave functions, 
see l53l. 
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Appendix D: Constructing tBLFQ basis: 
example 



An 



TABLE I. Spinor-polarization vector contraction for different 
helicity configurations of the incoming electron ("1"), outgo- 
ing photon ("2") and the outgoing electron ("3"). 



In this Appendix we illustrate the construction of the 
extended BLFQ basis, the diagonalization of the Hamil- 
tonian and the construction of the tBLFQ basis. For the 
sake of clarity we work in a highly truncated basis space 
for which numerical results are subject to large trunca- 
tion error, but this section is for illustration only. 



1. BLFQ state enumeration 

We take 6=me and L=27r MeV^^. In this ex- 
ample, our extended BLFQ basis consists of two 
segments, labeled by {^^=3/2, M,=l/2, 7V/=1} and 
{i^=5/2, Mj=l/2, 7V/=1}. In each segment we truncate 
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the transverse degrees of freedom at 7Vniax=2. 

Consider now which Fock states are present in our ba- 
sis. In each segment, the basis states have to meet the 
symmetry constraints (22 1 to (24), and the truncation 



constraint (28). Consequently, there are only two states 
in the K='i/2 segment; their quantum numbers are given 
in Table [TT} In the K=b/2 segment we have 3 basis states, 
see Table IIII) The total dimensionality of this extended 
BLFQ space, which is simply the sum of the two seg- 
ments, is 2+3=5. Note that we have assigned a number 



basis state no. 


Fock-sector 




n" 




A" 


k'< 




m'' 




1 




3/2 








1/2 










2 


1 n) 


1/2 








-1/2 


1 








1 



TABLE II. BLFQ basis states in the segment {if =3/2, 
M,=l/2, Nf^l}. 



basis state no. 


Fock-sector 




n" 


m" 


A" 


k'-' 




m"' 


A^ 


3 




5/2 








1/2 










4 


1 n) 


3/2 








-1/2 


1 








1 


5 


1 n) 


1/2 








-1/2 


2 








1 



TABLE III. BLFQ basis states in the segment of {K=5/2, 
M,=l/2, Nf=l}. 

to each of the basis states; this ordering is a matter of 
choice, and we include it so that the structures of the ex- 
act QED eigenstates which we will construct below, can 
be more easily related to their Fock components. 

Now that we have the basis states | a ) , we cal- 

I 



culate their matrix elements with the Hamiltonian 
( a' I- Pqed I ); using for example the expressions (C 
and ( C9 ) . The resulting QED Hamiltonian in our 5 



dimensional extended BLFQ basis is shown in Table |IV| 
The diagonal entries come from the free kinetic terms in 



(MeV) 


BLFQ basis state | a ) 


1 |e) 


2 167) 


3 le) 


4 67) 


5 |67) 


BLFQ 

basis 
state 
(a'i 


1 


^> 


0.3482 


-0.0119 











2 


1 67) 


-0.0119 


0.9139 











3 










0.2089 


-0.0024 


-0.0101 


4 


1 n) 








-0.0024 


0.3917 





5 










-0.0101 





0.8486 



TABLE IV. The Hamiltonian matrix -Pqed in the 
extended BLFQ basis consisting of two segments: 
{K=3/2,Mj=l/2,Nf=l} and {K=5/2, Mj=l/2, Nf=l}. 
See text for truncation parameters. 



the Hamiltonian. (In general, the kinetic terms give off- 
diagonal matrix elements, as can be seen from ( C8 ); such 



terms do not appear here only because of the small ba- 
sis space.) The off-diagonal matrix elements in Table IV 
come from the vertex interaction Wi. As expected, the 
Hamiltonian matrix exhibits a block-diagonal structure; 
no coupling exists between the K—3/2 segment (states 1 
and 2) and the K—5/2 segment (states 3, 4 and 5), due to 
the symmetries of QED. Normally, we would also include 
the instantaneous W3 terms. However, as noted above, 
we are only solving for QED mass eigenstates with the 
Hamiltonian accurate to order e in the present work. 



2. Diagonalization of -Pqed 



^QED eigenstate 

{P\ 


P^'(MeV) 


M^(MeV) 




BLFQ amplitudes {/3\a) 


K 


l|e) 


2|e7) 


3|6) 


4|e7) 


5 1 67) 


1 


0.3479 


0.5106 


3/2 


-0.9998 


-0.0210 











2 


0.9142 


1.0540 


3/2 


-0.0210 


0.9998 











3 


0.2087 


0.5105 


5/2 








0.9998 


0.0130 


0.0157 


4 


0.3917 


0.8474 


5/2 








-0.0130 


0.9999 


-0.0003 


5 


0.8488 


1.3640 


5/2 








-0.0157 


0.0001 


0.9999 



TABLE V. Eigenstates and eigenvalues of the QED Hamiltonian in the extended BLFQ basis comprising the two segments 
{A'=3/2,Mj = l/2,7V/ = l} and {K=5/2, Mj = l/2, Nf = l}. 



With the QED Hamiltonian matrix in the BLFQ basis prepared, we are ready to diagonalizc it. This can be done 
segment by segment because of the block-diagonal structure of the Hamiltonian. Doing so, we obtain eigenstates | (3 ) 
and eigenvalues . These are listed in Table. |vj The first column enumerates the eigenstates. The second and 
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third cohimns contain and the invariant mass for each eigenstate, which wiU be discussed shortly. The fourth 
column contains the segment specifier; it is enough to give just K in this case. In the 5th to 9th columns we list the 
overlaps of the eigenstates | /? ) with the BLFQ basis states | a ) in Tables [lT| and III i.e. this part of the table contains 
the coefficients in the expansion 

|/?)=^|a)(a|/3). (Dl) 

a 

Let us comment briefly on the physical interpretations of these QED eigenstates. In this example, the five QED 
eigenstates lie in two segments. The eigenstates numbered 1 and 2 are in the (total longitudinal momentum) K = 5/2 
segment, while eigenstates 3, 4 and 5 are in the K — 5/2 segment. In order to interpret these states it is useful to 
introduce the invariant mass, 

:= P+P- - P^P^ . (D2) 
We see that the eigenstates 1 and 3 have invariant masses close to the physical electron mass rrie. Reading off the 



coefficients (Dl), we see that these states are dominated by contributions from the single electron basis states. Thus 
we interpret them as the physical single electron states | Cphys ) with different longitudinal momenta. (The slight 
deviation in the invariant mass from me is due to our omission of mass corrections from counter-terms.) The small 
I e"f) components in their wavefunctions are generated by the QED vertex interaction (W^i), and describe the dressing 
of the bare fields by the photon cloud which, together, make up the physical electron [37H39] . We will see that these 
I 67) components play an important role in photon-radiation processes. They are also responsible for the electron's 
anomalous magnetic moment, see Ref. [JB] for more details. 

The eigenstates 2, 4 and 5 are excited states in their respective segments, with invariant masses considerably above 
the physical electron mass. Since they are dominated by the basis states in the | 67 ) sector, it is natural to interpret 
them as the electron-photon scattering states | e7scat ) • Their invariant masses = (Pg -t- P-y)^ are experimentally 
accessible through simultaneous measurements of the electron and photon four momenta. The | e7scat ) states receive 
small contributions from the single Fock electron sector ( | e ) ) due to the QED vertex interaction Wi . As we will 
see below, it is through such "minor" components that external fields are able to couple physical electron states to 
electron-photon scattering states or, in other words, cause photon emission. 



3. tBLFQ: time evolution 



The eigenstates we have constructed comprise the tBLFQ basis. We continue our example by calculating transitions 
between these eigenstates in nCs. The laser profile used in this example is 

eA~{x~) — 2aQmeCos{l_x~) — aQme[exp{il_x~) + exp{—il_x~)^ , (D3) 

and we take ao=l, which is at the edge of the nonperturbative intensity regime. Recalling that the frequency l_ can 
be written in terms of the wave number fcias as ^_ = ^fcias, we take fcias=l. The laser can therefore cause transitions 
between just the two segments of the tBLFQ basis prepared above. 

We need the matrix elements of V in the tBLFQ basis. One can first write down the matrix elements in the 
extended BLFQ basis, i.e. the set (a' |y| a), and then transform to the tBLFQ basis using 

(/3'I^I/3>=E(/5'I"')("'I^I")("I/^) • (D4) 

The resulting matrix of V in the tBLFQ basis, (/?' \V\ (3), is shown in Table. [Vl| The only allowed transitions are now 
between the two segments, because of the longitudinal momentum being added to the system by the background. The 
most probable transitions are those between the physical electron states | Cphys ) in the two segments, and between 
the electron-photon scattering states | e7scat ) in the two segments. These types of transitions describe acceleration, 
as the particle number is conserved, but the longitudinal momentum is changed by one unit. 

There are also transitions between the physical electron states and the electron-photon scattering states (| Cphys ) ^ 
|e7scat)), which describe the radiation process. We see that these transitions have much smaller amplitudes, since 
they link the "minor" Fock components in tBLFQ basis states {e.g., the \ ej) components in |ephys), and the 
I e) components in | e7scat )), which are suppressed by one factor of the electron charge e. Next we multiply by the 
required phase factor e^^'^'f'^ which transforms the matrix elements into those in the interaction picture, 

( \Viix+)\P ) = {l3'\V\P)x e^"^'^-V2 ^ 
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(MeV) 


Basis element | /? ) 




1 


2 


3 


4 


5 


Basis 
element 


1 (7^=3/2) 








-0.5109 


-0.0041 


0.0080 


2 {K=2,/2) 








-0.0041 


0.5110 


0.0002 


3 (^-=5/2) 


-0.5109 


-0.0041 











4 (^-=5/2) 


-0.0041 


0.5110 











5 (^"=5/2) 


0.0080 


0.0002 












TABLE VI. The matrix elements of the interaction term V in the tBLFQ basis. 



where w^'^ — Pg, — Pg and (see Table. |vj for the values of these energy eigenvalues). Due to this phase factor, 
the transition amplitudes oscillate in time with the period, ~l/w/3'/3, inversely proportional to the light-front energy 
difference between |/3') and |/3). Thus in cc"*"— >oo limit only the transitions which conserve energy can accumulate. 
The interaction picture matrix elements of Vj are given in Table. |VII[ We also need the initial state of the system. 



{P'\Vi{x+)\P) (McV) 


tBLFQ basis state | j3 ) 


1 


2 


3 


4 


5 


tBLFQ 
basis 
state 
1/3') 


1 (A"=3/2) 








-0.5109600^0"^ 


-0.0040e-oo22-+ 


0.0080e-0-250ia; + 


2 (^-=3/2) 








-0.0040eO-^'5^""^ 


0.5110e°-^""'^ 


0.0002e003^'"'^ 


3 (^-=5/2) 


-0.5109e-°°™^""^ 


-0.0040e-°-^^^""^ 











4 (^-=5/2) 


-0.0040e°°22-+ 


0.5110e-''-2''""'^ 











5 (^-=5/2) 


0.0080e°-2^°'""^ 


0.0002e-''"3='"'*' 












TABLE VII. The matrix elements of the interaction term Vi in the interaction picture, x'^ is in units of MeV 



I V'; a^"*" = 0)/, cf. Eq. In the nCs process the initial state is a physical electron. In our current tBLFQ basis there 
are two states corresponding to physical electrons, the two | Cphys ) states, see Tables O and 



III 



We choose our initial 



state to be that in the K = 3/2 segment. Now we are in a position to evolve the initial state forward in x+. To do 
so, we must identify the largest eigenvalue of Vi in order to determine our step size (see the discussion in S ect. IV C I. 
According to Table. VII the smallest value for l/ujp'p is about ^3MeV~^, and according to Tabic. VI the largest 
eigenvalue of Vj is about 0.511MeV which translates to 1/V7;max '-^2MeV~"'^. For this problem, we can safely choose 
the step size of Ja;+=0.1MeV~^ 





x+ (MeV-^) 





0.05 


0.1 


0.2 


1.0 


5.0 


10.0 


15.0 


20.0 


1 {K=Z/2) 


1.000000 


1.000000 


0.999347 


0.997391 


0.936139 


0.093892 


0.720013 


0.566033 


0.199441 


2 (^-=3/2) 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000001 


0.000010 


0.000005 


0.000005 


3 {K=h/2) 


0.000000 


0.000163 


0.000653 


0.002609 


0.063841 


0.905839 


0.279520 


0.432517 


0.798676 


4 (ii'=5/2) 


0.000000 


0.000000 


0.000000 


0.000000 


0.000004 


0.000060 


0.000035 


0.000005 


0.000015 


5 (^-=5/2) 


0.000000 


0.000000 


0.000000 


0.000001 


0.000016 


0.000209 


0.000422 


0.001440 


0.001862 




1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


1.000000 


0.999999 


(M^) (MeV) 


0.510635 


0.510635 


0.510635 


0.510635 


0.510643 


0.510745 


0.510985 


0.511826 


0.512153 



TABLE VIII. The evolution of I Gphys ) with K—3/2 in the laser field by Eq. ( |D3[ ). The norm and invariant mass are listed in 
the 6th and 7th row, respectively. 



Startin g from the initial state | ip; )/, i.e. a physical electron with K = 3/2, the evolved state | ^p; )/ is represented 
in Table. VIII at various times by the probabilities for being in various QED eigenstates | /3 ) . As time begins to 
evolve, the state first acquires an overlap with the | ephys ) state with K = 5/2 (the state with {f3,K} = {3,5/2}). 
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FIG. 8. Feynman diagram representation of the perturbative processes contributing to the matrix element of the background 
field vertex between exact QED eigenstates. 



Mathematically, this state is populated first due to the large matrix element in Vi coupling it to the initial state. 
Physically, this is acceleration; the electron is accelerated by the laser but does not, yet, have a significant probability 
for photon emission. 

At later times, the | e7scat ) states ({/?, if }={4, 5/2} and {5,5/2}) become populated. The overlap between these 
states and the initial state is smaller than between the physical electron state with the initial state. Hence, we observe 
a "domination" of acceleration over radiation, as the latter is suppressed by a factor of the coupling. Note that the 
I e7scat ) scattering state in the K = 3/2 segment, {{f3, K}={2,3/2}) eventually becomes populated, even though it 
is not directly coupled to the initial state. This population arises through the basis states in the K—5/2 segment, 
which are "decelerated" to the {/3, K} —{2, 3/2} state by the laser field. (This is due to the presence of the "negative" 
exponential in the chosen laser profile ( |D3[ ), which subtracts rather than adds longitudinal momentum.) 

The norm of the state \ip;x^) is conserved over the entire evolution due to the stability of the MSD2 scheme. It 
is an advantage of the tBLFQ approach that the wavefunction of the system is accessible at each time step, which 
allows one to monitor the real-time evolution of any observable 0{x'^), by evaluating 0{x~^)=i{il'', \Oj\ip;x^ )i. 
For example, taking Oi to be the invariant mass operator, we obtain the evolution of the average invariant mass of 
the system, see the 7th row of Table. VIII The increase in the invariant mass with time reflects the fact that the laser 
field pumps energy into the system and hence that photons are being created. 

This completes our example. Performing calculations in larger basis spaces follows the same procedure. In Sec. |VB| 
we present the numerical results for the nCs process in larger basis spaces, which allows for a more accurate description 
of the evolution of the system. 



Appendix E: Comparison with perturbation theory 

In order to check the BLFQ calculation, we calculate the matrix element ( e7scat \V\ Cphys ) in perturbation theory. As 
discussed in Sect.jVAj the perturbative approximation to this matrix element follows from standard time-independent 
perturbation theory, and is equal to 

( e7scat \V\ ephy. ) = ( 67 \Vq ^ V\e} + ( \V ^_ VqIc) . (EI) 



The Feynman diagram representation of the perturbative matrix elements on the right hand side of (El) is shown in 
Fig. [sj these are (respectively) s-channel and i-channel Compton scattering diagrams, where the laser field takes the 
place of the incoming photon. Using the machinery in Appendix [B] to write these perturbative overlaps in terms of 
the BLFQ basis, one eventually finds (with l^=2l_ from ([16]) and (5* the Kronecker delta). 



and 



(E2) 



^f;ee(e7) - ^frec 4V2L ] (2^)Vfc'^ 



(E3) 
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The final factors in eacfi of the above, describing spin and polarization contributions, can be read off from Table [T] 
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